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We derive an exact result for the averaged Feynman propagator and the corresponding density of 
states of an electron in two dimensions in a perpendicular homogeneous magnetic field and a Gaussian 
random potential with long-range spatial correlations described by a quadratic correlation function. 
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The problem of an electron moving in two dimensions 
(2D) in a strong homogeneous magnetic field B and a 
random potential has been studied particularly inten- 
sively since the discovery of the integer quantum Hall 
effect (IQHE) ||l[. Especially the density of states (DOS) 
as the simplest quantity that characterizes the quantum 
mechanical spectrum has been the subject of continu- 
ous interest P-p^. This was partly induced by the fact 
that it remains notoriously difficult to explain the un- 
expectedly large density of states between the Landau 
bands (LBs) detected experimentally |l^-[l^]. The theo- 
retical calculations, mostly based on models with short- 
range correlated randomness as, for instance, the Gaus- 
sian white noise potential, yielded in the high magnetic 
field limit well separated, disorder broadened LBs with 
a width proportional to ^/B and an exponentially small 
DOS in between the bands |l^Jl^. On the other hand, 
the experimental data indicate that there is a roughly 
constant background DOS which contains as many as 
10% to 20% of the levels |l|. 

It is very tempting to suspect that spatially long-range 
correlations of the disorder, due to inhomogeneities, for 
instance, are responsible for this background DOS. Such 
a model was studied by using path integrals in connec- 
tion with a cumulant expansion . This reproduced the 
experimental result. However, the infiuence of the cu- 
mulant approximation which is uncontrolled in the pres- 
ence of a magnetic field was unclear. Subsequent ap- 
proximative treatments of the DOS in the presence of 
long-range correlated disorder, by using functional inte- 
grals and without the cumulant approximation, did not 
improve the situation Exact results for long-range 
correlated disorder exist only for the band tails |j . 

In this paper, we consider a model that has already 
been treated before without magnetic field |2^. Also in 
the region of the quantum Hall effect, this model allows 
for the treatment of the infiuence of long-range correla- 
tions in the randomness. The essential assumption is a 
quadratic approximation of the correlation function. The 
range of validity of this model assumption is discussed 
below. The Feynman propagator can be calculated ex- 



actly, and the DOS evaluated numerically, but without 
employing further approximations. 

The comparison with recent experimental results pl| ] 
shows that long-range correlations can qualitatively ex- 
plain some aspects of the experimental data. However, 
the quantitative behavior, as a function of the magnetic 
field, cannot be understood. We conclude, that the large 
DOS between the LBs seen in experiment, is very likely 
related to interaction between the electrons. 

We consider a spinless electron moving in 2D under the 
infiuence of a homogeneous magnetic field and a random 
impurity potential V{r). The Hamiltonian is 
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H = - [iV + A{r)Y + V{r) 



(1) 



For the vector potential, we assume the symmetric gauge 
A{r) — e-r/2, where e denotes the antisymmetric (2 x 2) 
matrix, so that = —1. The magnetic field B defines 
the intrinsic energy, length, and time scales. Energies and 
lengths are measured here in units of the cyclotron energy 
hujc = TieB/m and the magnetic length = &B, 
respectively (e elementary charge and m effective mass). 
The unit of time is then the inverse of Wc- The probability 
distribution of the impurity potential Vir) is assumed to 
be Gaussian with zero mean (V^(t)) = and correlation 
function 



C(r,r') 



iy{r)V{r')) = 



1 



1 



(2) 



Angular brackets denote the ensemble average with re- 
spect to the probability distribution of the random vari- 
able V . The parameters r\ and L determine the strength 
and the range of the impurity potential, respectively. 

The above polynomial correlation function Eq. (y) is 
the first order expansion of a Gaussian 

CG(r,r') = 7y2e-(^-^')Vi^ (3) 

The function C(r,r'), Eq. (^, has been used previously 
pof to describe approximately the DOS of an electron 
in 2D in a random impurity potential with the Gaussian 
correlation function, Eq. (||). 
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Here, we derive the exact result for the impurity aver- 
aged Feynman propagator of the model Eq. (||), 

r{t)=r 

{K{r,t;ro))^{ J V[r]exp 
r(o)=ro 



J I?[r] exp i J dT£{r, r; r) 



(4) 

where I?[r] is the Feynman measure on the set of trajec- 
tories with boundary conditions r(0) = ro and r{t) = r. 
The time evolution begins at t = 0, which is omitted 
in K{r, t; ro) for simplicity. Equation is a Feynman 
path integral and the Langrangian 



L{r,r;T) = -r^ 



-r er 



V{r) 



(5) 



represents the usual expression for a charged particle in 
2D in a magnetic field and a potential. 

Since V{r) enters only linearly in the exponent of 
Eq. (0), the impurity average is easily performed, 



r(t)=r 

{K{r,t;ro)) ^ J V[r]exp 
r(o)=ro 



dr- 

2 L 



X exp 



dr / dr' C{r{T),r{T')) } . (6) 
Jq 



The last factor in Eq. (|6|) results from the impurity av- 
erage. It is purely real and non-local in time. This con- 
sequence of integrating out degrees of freedom is well 
known ||2^. It is usually the end of an exact calculation. 
In our model, Eq. (|^), however, we have the simplifying 
feature that C {r,r') is quadratic in r — r'. Therefore, 
the path integral remains quadratic and can be evaluated 
exactly by using standard techniques - in spite of the 
non-locality in time. Details are given in p^ ]. 
Denoting {K{r,t;ro)) = KB{r,t;ro), we obtain 



KB{r,t;ro) 
1 



(iV4)(f^2 _ 1/4) 



2Trit sin [{Q + 1/2) {t/2)] ■ sin [{n - 1/2) (t/2)] 
1 



and 



/(r,ro;0 



X exp 



n 

2 sin fit 



f]H' + ifir,ro;t) 



(7) 



{r " ro) 



■ cos fit 



-2r^e^'*r + ^ ^ cos -t 

° 2 2 



{r-rof {nsin^t- isinffi)' 



+ 



2172 



2rn er 



cos it — cos fit 



sin —t 
2 



■ sin Qt 



(8) 



The complex frequency = ^^/l + 8it{ri/L)^ depends 
on the time t of propagation. Configurational averaging 
leads to a path integral with a non-local action. There- 
fore, the above is not a conventional Feynman propaga- 
tor. It does not have the group property, 

KB{r,t;ro)^ J d\' KB{r,t - t' ■,r') Keir' ,t';ro), (9) 

which would be valid for a unitary time evolution opera- 
tor without the configurational average. 

This has serious consequences for the analytical struc- 
ture of the diagonal part 



KB{t) = KB{r,t-r) = 



2TTit 



(<V4) _ 1/4) 



. [{n + 1/2) {t/2)] sin [(17 - 1/2) {t/2)] 



(10) 



of our propagator that determines the DOS. By investi- 
gating the analytical continuation of KB{t) into the com- 
plex t-plane, we find not only a single pole for t = on 
the real time axis, but also poles in the lower and up- 
per complex plane. In the limit L oo, the poles in 
the lower plane approach the real time axis whereas the 
poles in the upper plane disappear towards infinity. 

Evaluating the DOS for such a propagator, we have to 
take two conditions into account: first, causality should 
be fulfilled and second, the definition of the DOS for this 
model should reproduce the limit L = oo of a model with 
constant spatial correlations. The DOS can be calculated 
exactly in this limit [^,^. We choose the integration 
contour in the complex t-plane in such a way that coming 
from minus infinity all poles, t„, in the lower half plane 
are passed anti-clockwise. This choice is unique, and the 
limit L — > oo is continuous. Deforming this contour back 
to the real axis gives the DOS 



D{E) 



An 



oo 

i-P j dtKB{t)e' 



Et 



+^ Res,=,„ [ifB(t)e^^'] . (11) 

lm(t„)<0 

Here, V dt . . . denotes Cauchy's principal value of 
the integral and Rest=t„ [. . .] the residue at t = t„. 

Equation ([ill) is our central result. Figure 1 shows the 
result of the integration for L — 6 and rj — 0.2, 0.3 and 
0.4 (in the units defined above), normalized to the DOS 
at i? = 0, Dq. Due to the impurity potential, the Lan- 
dau peaks are broadened as expected. The broadening 
decreases for increasing energy which indicates that dis- 
order becomes less effective. For large L, Eq. ( [ll| ) can 
be shown to yield the same as obtained for the model 
with a constant correlation function, a superposition of 
Gaussians centered at the Landau levels (LLs). How- 
ever, Eq. (O) becomes negative if E exceeds a critical 
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value. This is a consequence of the unphysical shape of 
the polynomial correlation function Eq. (|^), as we will 
now see. 

While the correlation function Eq. (||) is positive def- 
inite, its truncated expansion Eq. (^) is not. Neverthe- 
less, using Eq. (^) in calculating physical quantities may 
give the same result as using Eq. (^) and expanding af- 
terwards in L~^. It is instructive to take the limit in 
Eq. (|l^) in which one keeps only the lowest LL in Eq. (|^) . 
Formally, this is achieved by measuring the energy from 
the center of the lowest LL, scaling E — fiLOc{E — 1/2) 
and f) = fiLiJc'ri, and taking lOc — >■ oo. We find for the DOS 
of the lowest LL with the random potential described by 
Eq. (i), 



Dlll{E) 



1 

2^ 



dt 
2^' 



2{fjt/Lf cos Et 
l-exp{-2{fit/Ly 



(12) 



The limit L — > oo is correctly reproduced. For L > \/2, 
Dlll{E) is positive for all energies. Now, the follow- 
ing argument suggests that including more and more LLs 
yields an increasingly higher bound for L until it diverges 
in the limit of all LLs taken into account. The classical 
cyclotron radius is proportional to y/E. Thus, the cor- 
relation function Eq. (|) is positive and the DOS well 
defined only for energies smaller than i^max ~ (^/2)^, 
independent of the strength of the disorder rj. Beyond 
this critical energy, the classical cyclotron radius, which 
is a characteristic length for the range of the electron, is 
of the order of L and larger. Then, the expansion lead- 
ing from Eq. (||) to Eq. (||) breaks down and, for such 
distances, the correlation function Eq. (||) becomes neg- 
ative. Such a bound -Emax, and its independence of rj, 
can indeed be seen in Fig. 1. 

Our model is the first order approximation of an im- 
purity potential with Gaussian correlations, Eq. (|^), 
which is believed to be typical for such materials as 
GaAs/AlGaAs heterostructures In order to obtain 
information on the validity of our model, Eq. (H), we 
compare the DOS with previous results for the Gaussian 
model. 

By neglecting the mixing of states between differ- 
ent LLs and assuming a Gaussian shape of each level. 



Broderix et al. 



found for the width of the 



Pr 



V L4 



LB 



(13) 



with Pn{x) the n^^ Legendre polynomial. Figure 2 shows 
both results, ours and that of |2^], for L = 10 and rj = 0.4. 
The agreement is quantitative. For increasing L, the 
agreement becomes even better. This demonstrates in- 
deed that for large values, the length parameter L has the 
meaning of the correlation length of a Gaussian correla- 
tion function. Also, our results support the conjecture 



that in the quantum Hall regime, LL mixing is negligible 
and the DOS is a superposition of Landau bands with 
Gaussian shapes. This is valid for even a broader energy 
range when L is increased. 

Armed with the exact result for the DOS of a model 
with spatially correlated randomness, we now compare 
with experiments, cf. Fig 3. By plotting the DOS at 
the cyclotron energy - between the lowest and the sec- 
ond lowest LL - against the magnetic field strength, we 
find that the various experiments on GaAs/AlGaAs het- 
erostructures show the same qualitative behavior. How- 
ever, this appears to be quite different from the predic- 
tion based on the model of non-interacting electrons with 
spatially correlated disorder. We find a strong decrease 
of the DOS when increasing the strength of the mag- 
netic field, much stronger than is observed experimen- 
tally. Since a spatially correlated random potential en- 
ergy represents the most general model of disorder, this 
leads to the question, whether or not the interaction be- 
tween the electrons has to be taken into account, even 
for a theory of the DOS in the region of the IQHE which 
is commonly believed to be explainable without electron- 
electron interactions. More quantitative comparisons be- 
tween experiment and theory are necessary, in order to 
clarify this point. 

If our above suspicion proves to be true, the proper 
understanding of the IQHE would very probably need, 
in addition to disorder, electron-electron interaction as 
an important ingredient, as is already the case for the 
fractional quantum Hall effect. 
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FIG. 1. DOS, in units of the DOS at B = 0, for the eight 
lowest LLs for L = 6 and different strengths of disorder: thin 
solid line rj = 0.2, dotted line rj — 0.3, thick solid line rj — 0.4. 



FIG. 2. Thick solid line: DOS of the present model; thin 
solid line: approximation for a model with Gaussian correla- 
tions for L = 10 and 77 = 0.4. 



FIG. 3. Normalized DOS at the cyclotron energy 
versus the magnetic field strength. Solid line: experi- 
ment of [Q; solid line with bullets: recent experiment 
of Results for the present model: ri = 0.67 meV; 

L = 1089-10-*cm(A), = 1.26 meV, L = 787 • 10~**cm (Q), 
rj = 1.86 meV, L = 485 ■ 10"** cm (O). 
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